Isospectrality of conventional and new extended 
potentials, second-order supersymmetry and role of VT 

symmetry 



B. Bagchi a , C. Quesne R. Roychoudhury c 

a Department of Applied Mathematics, University of Calcutta, 
92 Acharya Prafulla Chandra Road, Kolkata 700 009, India 
b Physique Nucleaire Theorique et Physique Mathematique, Universite Libre de Bruxelles, 
Campus de la Plaine CP229, Boulevard du Triomphe, B-1050 Brussels, Belgium 
c Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700 035, India 

Abstract 

We develop a systematic approach to construct novel completely solvable rational 
potentials. Second-order supersymmetric quantum mechanics dictates the latter to 
be isospectral to some well-studied quantum systems. VT symmetry may facilitate 
reconciling our approach to the requirement that the rationally-extended potentials 
be singularity free. Some examples are shown. 
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1 Introduction 



The search for exactly solvable potentials has been a subject of intense research since the 
advent of quantum mechanics because they may provide a conceptual understanding of 
some physical phenomena, as well as a testing ground for some approximation schemes. 
This has continued to be so after the generalization of quantum theory for the case of 
complex PT-symmetric [1] or pseudo-Hermitian [2] Hamiltonians. 

Over the years supersymmetric quantum mechanics (SUSYQM) has emerged as one of 
the most insightful tools towards construction of Hamiltonians with a prescribed spectrum 
starting from some given exactly solvable form [3, 4, 5, 6]. The key technique such as 
the factorization method (sometimes couched in the language of intertwining relationships 
amongst operators inducing factorization) has enabled one to uncover many useful proper- 
ties of quantum mechanics. SUSY systems are appropriately described by a superpotential 
descendant of the given Schrodinger potential, which in turn furnishes the so-called SUSY 
partner Hamiltonians. In one-dimensional systems, strict isospectrality may be present in 
the first-order systems (broken SUSY) when the factorization energy E is smaller than the 
ground-state energy of the starting Hamiltonian [7, 8]. In such a context, the potential 
and eigenfunctions of the partner Hamiltonian are known in terms of analytic expressions, 
which, however, may be a good deal more complicated than the corresponding ones of the 
initial Hamiltonian [9] . Interestingly one can look for extensions to higher-order SUSY the- 
ories by resorting to higher-derivative versions of the factorization operators. These offer 
the scope of relating to non-trivial quantum systems as has been found in a number of 
recent works [10, 11, 12, 13]. 

Against this background, the recent introduction of two new classes of exceptional or- 
thogonal polynomials [14] and their occurrence in the bound-state wavefunctions of two 
novel rational potentials isospectral to some well-known conventional ones [15] have led us 
to re-examine the construction of such pairs of partner potentials. These examples have 
suggested us an alternative approach to the usual one, which consists in searching for a 
reduction of the initial Schrodinger equation general solution to some elementary function 
[16]. 
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Taking advantage of such an experience, we plan to enquire here how flexible is the solu- 
tion of the superpotential that forms a given SUSY quantum model. Employing the factor- 
ization method, which was generalized by Gendenshtein [17] in terms of the reparametriza- 
tion of the potential ensuring thereby an algebraic construction of the latter, we will gen- 
erate new pairs of SUSY Hamiltonians by effecting a reparametrization of the starting 
Schrodinger potential. We will show how the family of such pairs can be substantially en- 
larged should we redefine the existing superpotential in terms of new couplings and allow 
for the presence of some rational function in it. The latter plays an interesting dual role: it 
not only allows coincidence of one of the new partner Hamiltonians with the starting one 
up to some reparametrization of couplings, but also forces the other partner to emerge as a 
superposition of the initial exactly solvable potential (with its parameters left undisturbed) 
and some additional rational terms, thus representing an extended family of exactly solv- 
able rational potentials. Such an extended potential has its spectrum unaltered vis-a-vis the 
exactly solvable one from which it is obtained, while the corresponding bound-state wave- 
functions turn out to be expressible in terms of rather simple polynomials. It is important 
to stress that a second-order supersymmetric (SSUSY) transformation lies at the very root 
of the observed isospectrality. Furthermore, our findings show that for some potentials VT 
symmetry [1] enters the picture in a natural way to render the derived potential singularity 
free. 

2 Generalized Poschl- Teller potential 

We begin with the concrete example of the generalized Poschl- Teller (GPT) potential as 
our first model to address the issue of constructing an extended class of rational potentials. 

2.1 Extended family of rational potentials 

The GPT potential reads in standard notations [3] 

V AB (x) = [B 2 + A(A+l)]csch 2 x - B(2A+l)cschxcothx (1) 
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and is denned on the half-line < x < oo. It is repulsive at the origin for B > A + 1 > 1 
but supports a finite number of bound states, whose energies are given by 



E \f) = -{A-v)\ z/ = 0,l,...,z/ max (A-l<u max <A). (2) 

The associated wavef unctions, vanishing at the origin and decaying exponentially at infinity, 
can be expressed in terms of Jacobi polynomials as 

^ B \x) =AC (AB) (cosh:r- l)^ B ~ A \coshx + l)-^ B+A ^p!) B ~ A ~ h ~ B ~ A ~ h \coshx), (3) 

where 



kt(a,b) = 2 a( u\(2A-2u)T(B + A-u + l) 
\r(B-A + v+±)r(2A-v + l) 



1/2 



denotes a normalization constant 1 . 

In the standard SUSYQM approach to such a potential, which incidentally is also shape 
invariant, deletion of the ground-state energy yields a partner potential Va-i,b(x) of 
the same kind by using a superpotential having the form A coth x — B csch x [3] . 

Motivated by the previous results on the non-uniqueness of factorization schemes which 
also have relevance to the construction of rational potentials [7, 16, 18, 19, 20], let us 
now ask the following question: can we modify the above superpotential for the GPT 
potential in such a way that the starting potential (x) = W 2 — W + E (E denoting the 
factorization energy) continues to belong to the GPT family up to some reparametrization 
of coefficients, while its partner V^\x) = W 2 + W' + E differs from Va,b(x) in the presence 
of some additional rational terms? Taking cue from the observation that z = cosh x is the 
basic variable appearing in the Jacobi polynomial that controls (3), we propose for such a 
superpotential the following form 

TTr , , , , i sinhx . . 

W(x) = a coth x + ocschx , (4) 

v ' coshx + c v ' 

where the constants a, b, c are to be determined from suitable inputs. 



1 Such a picture is also valid whenever the parameters vary in the range A + 1 > B > A + \ > \ 
or A + I; > B > A > 0, where the potential is weakly attractive at the origin. Among the two square 
integrable solutions of the Schrodinger equation corresponding to a given , it is customary to consider 
the function (3) as the physical one because it is the most regular. 
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Straightforward calculations using (4) lead to 
\/(+) ( x ) = [a(a + l) + b 2 ] csch 2 x + (2a + 1)6 csch x coth x + + ^ + ( ° ~ + E (5) 



and 

V (_) (a:) = \a(a - 1) + 6 2 1 csch 2 x + (2a - 1)6 csch a; coth a; + ^ — — 
L coshx + c 

2(c 2 - 1) | | £ (6) 

(coshx + c) 2 

Elimination of the term proportional to (coshx+c) -1 from V^ + \x) can be achieved provided 
we choose c = 26/ (2a— 1). Furthermore, the first two terms of V^~\x) coincide with Va,b(x), 
given in (1), if the two conditions 

a(a-l) + 6 2 = J B 2 + A(A + l), (2a - 1)6 = -B(2A + 1) (7) 

are fulfilled. Finally, the vanishing of the constants terms in (5) and (6) imposes that the 
factorization energy be given by E = — (a — l) 2 . 

The two restrictions in (7) lead to four possible solutions for a and 6. The first two are 
given by a = ± (5 ± |) , 6 = =F (A + |) , where we take all upper or all lower signs, while 
the two remaining ones can be obtained from them by permuting A+ 1 with B. As a result, 
c becomes c = — (2 A + 1)/(2B) or c = — 2B/(2A + 1), respectively. Only for the former 
c value can we assert that the extended potential V^~\x) will be free from poles on the 
half-line. We therefore restrict ourselves to such a value and arrive at the following results: 

V {+ \x) =Va,b±i(x) = [(fi±l) 2 + A(A + l)]csch 2 x-( J B±l)(2A + l)cscha;cothx, (8) 

2(2,4+1) 2\4B 2 - (2A + 1) 2 1 



with 



W(x) = ± [B ± - cothrr t U + - csch a; - 



2/ \ 2/ 2 J Bcoshx -2A-V 

E = -[B T 



V 2 



(10) 

We have therefore found two distinct solutions to the above-mentioned question: going 
from a conventional GPT potential to an extended potential VA,B,ext(x) can be achieved in 



a first-order SUSYQM setting provided we start either from Va,b+i(x) or from Va,b-i(x) 
(but in no way if we start from a potential Va,b{ x ) with the same parameters). It is also 
worthwhile to note that although the GPT potential Va,b(%) is not defined for B — A + |, 
the same is not true for VA,B,ext(x), which can easily be shown to be equivalent to the well- 
behaved conventional GPT potential V A _ 1A+ z{x), with a correspondence between their 
respective bound-state energies and wavef unctions. 

The facts that the GPT energy eigenvalues (2) are independent of B (provided this 
parameter satisfies the condition B > A + 1) and that the factorization energy E, given 
in (10), is smaller than make it plausible that the partner potential V <y ~\x) in (9) 
be strictly isospectral to any choice for V^ + \x) in (8) (as well as to Va,b(x) in (1)). This 
is actually confirmed by determining the factorization function (f>(x), which is such that 
W = —(p' /(p. The result reads 

4>{x) oc (cosh a; - l) T ^ B ' A± ^) (coshx + 1) T ^ B+A± ^) (2B coshx - 2A - 1), 

where the upper (resp. lower) signs correspond to Va,b+i( x ) (resp. Va,b-i(x)) in (8). We 
note that for the former choice, <f>(x) decreases exponentially for x — > oo, but increases 
as x~( B ~ A ^ for x — > 0, while for the latter choice it vanishes as x B ~ A ~ l for x — > 0, but 
increases exponentially for x — > oo, thus showing in both cases that neither nor 0" 1 is 
normalizable, as is required by the general theory [8]. 

2.2 Determination of bound- state wavef unctions 

The bound-state wavefunctions ipi ■* (x) = ifj( A ' B ' cxt ^ ( x ) of the extended GPT potential can 
be determined by using the standard intertwining relations H^A* = AH^~\ AH^ = 
H^A [3, 4, 5, 6], satisfied by the partner Hamiltonians 

#(+) = A* A = — f-, + V {+ \x) - E, = ii f = — + V { -\x) - E. 

dx z dx z 

Here A* = — + W(x) and A = ^ + W(x). Hence choosing, for instance, Va,b+i(x) in 
(8), we can write 

^-\x) = ^-A^\x), 4 + \x) = ri A > B+1 \x), e u = E^ - E, (11) 
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for any v = 0, 1, . . . , z/ max . 

On using Eq. (10), as well as Eq. (3) with B + 1 substituted for B, we can rewrite (11) 

as 



where we have defined A/'i + ' ) = Af!) A ' B+1 \ z = z(x) = cosh a;, a — B — A—^, (3 — —B — A — \ 
and 

= [p + a -(P- a)z\ {{z - 1)^ + a + lj + {(5 - a)(z - 1). 

The action of the first-order differential operator oi a '^ on the Jacobi polynomial 
p(a+i,p i)^ can ^ j n f errec j f rom some known differential and recursion relations for the 
latter [21]. The result can be written as 

(*,p )p (a + i,p-i) {z) = _ 2{a _p ){u + a )P^\z), (13) 
where P^\z) is a [y + l)th-degree polynomial, defined by 

ffffw = -k* - t)p^(') + 1,: fl(2) - < 14 > 

2 a + p + 2z/ p — a 

in terms of some Jacobi ones. 

At this stage, SUSYQM can be profitably used to construct orthonormal bound-state 
wavefunctions ipi \x) in terms of P^\z). Inserting (13) in (12) and taking all previous 
definitions into account indeed leads to 

=JVJ">! 2 LsJ- 2.4-1 P i» "W*), < 15 > 



where 



= -45 - ^ A/" (+) 

Ju ^\ B + A-u-\) Ju 

= _ 2 a + 2 B ( * (2A-2v) {B + A-v+\)T(B + A-v-\) X 1/2 
(B-A + v+\)T(B-A + v-\)T{2A-v+l) 
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2.3 Isospectrality and SSUSY 

Let us now focus on the isospectrality issue of the conventional and extended potentials 
with the same parameters A, B. Since we have been able to go from Va,b±i{%) to VA,B,ext( x ) 
by a first-order SUSYQM transformation, we may try to convert Va,b{x) into Va,b±i(x) 
in a preliminary step. As is easily seen, such a transformation can be carried out in a 
first-order SUSYQM setting by assuming 

(+) = &B = — f- + V i+ \x) - E, H = B& = — f= + V { ~\x) - E, 

ax 2 ax 2 

V^\x) = Va,b( x ) = W 2 -W' + E, V ( -\x) = V a ,b±i(x) = W 2 + W' + E, (16) 

& = -4- + W(x), B = ^- + W(x), 
ax ax 

with a superpotential and a factorization energy given by 

W(x) = T{B±7))cothx±(A+±)cschx and E = -{B±\) 2 , 
respectively. As before, E < E^ and the corresponding factorization function 

4>(x) oc (cosh 3 :-l) ± K^ ± ^l)( cos ha; + 1)^(^+^+1) 

is such that neither nor _1 is normalizable, thereby ensuring the (strict) isospectrality 
of this step too. 

Since from (8) and (16), it follows that V^~\x) = V^(x), the two first-order SUSY 
systems ^H^ + \H^~^ and (H( + \ if( - )) can be glued together so as to get a reducible 
SSUSY one /i (2) ) [10, 11, 12, 13]. In such a framework, the two Hamiltonians and 
with corresponding potentials V^(x) and V ly2 \x) ) respectively, intertwine with some 
second-order differential operators 

Ai = ~ 2p{x) Tx + q{x) > A = &* + 2p{x) l + 2p ' {x) + q{xl (17) 

as Ah^ = h^A and A^h^ = h^A\ so that the functions p(x), q(x) and the potentials 
V^ 1,2 \x) are constrained by the relations 

q{x) = -p +p - — + 7T- 



2p \2pJ 16p 2 ' 
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16p 2 ' 



where c is some integration constant. 

The relation between both approaches follows from the equations 

&(!) = #(+) + £, = H^- |, A = &A\ A = AB, (18) 

where and /i( 2 ) are both partners of some intermediate Hamiltonian 

h = + - = Jf « - - 
2 2 

and the constant c is related to the two factorization energies through 

c = E-E = t2B. 

On comparing Eq. (17) with the two last relations in (18), we finally obtain 

1 /„, ~\ fisinhx 



P {x) = -(W + W) 



2B coshx-2A- 1' 

thereby completing the determination of the SSUSY transformation. It should be stressed 
here that since we have two possibilities for V^\x) = V( + \x) (see Eqs. (8) and (16)), there 
are two different paths for going from Va,b{x) to VA,B,ext(x) or, in other words, the SSUSY 
transformation admits two different decompositions and therefore two different intermediate 
Hamiltonians h. 

3 Generalizations 

The construction of extended potentials carried out in Sec. 2 in the GPT case can be 
generalized to some other exactly solvable potentials. 

3.1 Scarf I potential 

For the Scarf I potential 

V a ,b(x) = [A(A-l)+B 2 ]sec 2 x-B(2A-l)secxtcinx, -^-<x<^, < B < A-l, 

(19) 

for instance, the bound-state wavefunctions being expressible in terms of Jacobi polynomials 
in the variable z = sin a; [3], instead of (4) we would take W(x) = a tan a; + b sec x — 
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cosx(sinx + c)' 1 and obtain c = — (2 A — 1)/(2B) for an extended potential without pole. 
As a consequence, we would get V^(x) = Va,b±i(x) while its partner V^~\x) would 
coincide with the extended Scarf I potential obtained in [15]. The isospectrality of the 
latter with Va,b(x), given in (19), would then be explained by a SSUSY transformation 
entirely similar to that derived for the GPT potential. 

3.2 Scarf II potential and role of PT symmetry 

We may also inquire into what would happen for the corresponding hyperbolic Scarf II 
potential 

V a ,b(x) = [B 2 -A(A+l)}sech 2 x+B(2A+l)sechxtanhx, -oo < x < oo, A> 0, (20) 

for which the corresponding variable z is given by z = i sinh x [3] and we would therefore 
assume W(x) = atanhrr + frsecha; — i cosh x (i sinh x + c) _1 . Due to the fact that sinhx 
takes all real values whenever x runs over the real line, neither of the two solutions for c, 
namely c = — i{2A + 1)/{2B) and c = 2\Bj{2A + 1), would lead to an extended potential 
free from pole. 

It is interesting to note that this problem can be easily coped with by going to the 
PT-symmetric Scarf II potential, obtained by replacing B by LB in (20) [22, 23, 24, 25, 26]. 
Depending on the value assumed for c, this would result in two distinct "PT-symmetric 
extended potentials, the first one being given by 

V A ,iB, ex t(x) = -[B 2 + A(A + 1)] sech 2 a; + iB(2A + 1) sechrrtanha; 4 



2A + 1- 2\B sinh x 
2[(2A + l) 2 -45 2 ] 



(2A+ 1 - 2i J Bsinhx) 2 

and the second one following from it by permutation of A + \ with B. To each of them, 
there would correspond a different SSUSY construct. Such a doubling of extensions hints 
at a connection with the existence of two series of real energy levels, —(A — u) 2 , v < A, 
and — (5 — | — u) 2 , v < B — |, which characterizes the PT-symmetric Scarf II potential 
whenever A and B — \ are both positive [22]. 
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4 Concluding remarks 

In conclusion, we have devised a SSUSY approach to isospectral conventional and rationally- 
extended exactly solvable potentials, which, at the same time as providing us with a con- 
structive method for determining the latter, gives us a clue to calculating its polynomial 
solutions through the use of first-order differential operators. Such a procedure has been 
illustrated with the detailed example of the GPT potential and its application to the Scarf 
I potential of [15] has also been sketched. Furthermore, we have shown in the Scarf II 
case how VT symmetry may greatly help us in reconciling our approach with the require- 
ment that the rationally-extended potential be singularity free in the domain where the 
conventional potential is defined. 

In this work, for simplicity's sake, we have restricted ourselves to some shape invariant 
potentials, whose eigenvalues are independent of one parameter, and we have chosen to 
change this independent parameter to reparametrize the potentials. It should be stressed 
that, as we plan to show elsewhere, our method is by no way limited to such cases. 

A further important point, which has not been dealt with here for lack of space but 
will also be considered in a future work, is the fact that the polynomials appearing in the 
rationally-extended potential bound-state wavefunctions (see Eqs. (14) and (15)) belong 
to a class of exceptional orthogonal polynomials, which has been the topic of some recent 
mathematical study [14]. Our SSUSY approach provides us with a convenient way of 
constructing and generalizing such polynomials. 
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